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Abstract. 

We present the results of an investigation of the sequence of events from the Sun to 
the Earth that ultimately led to the 88 major geomagnetic storms (defined by minimum 
Dst < 100 nT) that occurred during 1996 - 2005. The results are achieved through 

cooperative efforts that originated at the Living with a Star (LWS) Coordinated Data- 
Analysis Workshop (CDAW) held at George Mason University in March 2005. Based on 
careful examination of the complete array of solar and in-situ solar wind observations, 
we have identified and characterized, for each major geomagnetic storm, the overall solar- 
interplanetary (solar-IP) source type, the time, velocity and angular width of the source 
coronal mass ejection (CME), the type and heliographic location of the solar source re- 
gion, the structure of the transient solar wind flow with the storm-driving component 
specified, the arrival time of shock/disturbance, and the start and ending times of the 
corresponding IP CME (ICME). The storm-driving component, which possesses a pro- 
longed and enhanced southward magnetic field (B s ). may be an ICME, the sheath of 
shocked plasma (SH) upstream of an ICME, a corotating interaction region (CIR), or 
a combination of these structures. We classify the Solar-IP sources into three broad types: 

(1) S-type, in which the storm is associated with a single ICME and a single CME at 
the Sun; (2) M-type, in which the storm is associated with a complex solar wind flow 
produced by multiple interacting ICMEs arising from multiple halo CMEs launched from 
the Sun in a short period: (3) C-type, in which the storm is associated with a CIR formed 
at the leading edge of a high speed stream originating from a solar coronal hole (CH). 

For the 88 major storms, the S-type, M-type and C-type events number 53 (60%), 24 
(27%) and 11 (13%), respectively. For the 85 events for which the surface source regions 
could be investigated, 54 (63%) of the storms originated in solar active regions, 10 (12%) 
in quiet Sun regions associated with quiescent filaments or filament channels, and 11 (13%) 
were associated with coronal holes. Remarkably, 10 (12%) CME-driven events showed 
no sign of eruptive features on the surface (e.g., no flare, no coronal dimming, and no 
loop arcade, etc), even though all the available solar observation in a suitable time pe- 
riod were carefully examined. Thus, while it is generally true that a major geomagnetic 
storm is more likely to be driven by a front-side fast halo CME associated with a ma- 
jor flare, our study indicates a broad distribution of source properties. The implications 
of the results for space weather forecasting are briefly discussed. 


1. Introduction 

A NASA Living With a Star (LWS) Coordinated Data 

Analysis Workshop (CDAW) was held at George Mason Uni- 

versity, Fairfax, VA, in March 2005. The workshop focused 

on the major geomagnetic storms of solar cycle 23, specifi- 

cally the 88 events from 1996 (corresponding to the start of 

observations from the SOHO spacecraft) to the end of 2005 

having minimum Dst (disturbance storm time index) < - 

100 nT. Four working groups were established to address (1) 

the solar and interplanetary (IP) sources of these storms, (2) 

storm mechanisms, (3) the associated ionospheric storms, 
and (4) storm predictions. Here, we summarize the efforts 
of Working Group 1 to identify the sequence of Sun- to- Earth 
activities for all 88 storms. The aim was to produce as com- 
prehensive a fist of solar-IP sources as possible by combining 
a wide variety of data sets and exploiting the different ar- 
eas of expertize of the group members. The purpose of this 
paper is to describe the identification methods and present 
the identification results, which we hope will serve as a ba- 


sis for further in-depth studies of these important Sun-Earth 
connection events. 

It is now well established that a geomagnetic storm is the 
consequence of a chain of causative events originating from 
the Sun’s corona and ultimately evolving into a geo-effective 
solar wind flow in near- Earth space [e.g., Brueckner et al. 
1998; Webb et al. 2001; Berdichevsky et al. 2002; Zhang 
et al. 2003; Gopalswamy et al , 2005]. Such geo-effective 
solar wind flows fall into two broad types, depending on 
their origins. One type is associated with an IP coronal 
mass ejection (ICME, also known as ejecta), the interplan- 
etary counterparts of CMEs at the Sun, and includes the 
disturbed shock sheath (SH) region upstream of the ICME 
(which may have a shock at the leading edge) and the ICME 
itself. The second type is associated with fast solar wind 
emanating from solar coronal holes, in particular with the 
corotating interaction regions (CIRs) that form at the lead- 
ing edges of such streams as they interact with the pre- 
ceding slower ambient solar wind. Previous studies have 
found that major/intense geomagnetic storms (e.g., Dst < 
-100 nT, or Kp > 7-) are mainly caused by ICMEs, while 
moderate and minor storms can be caused by both ICMEs 
and CIRs [Gosling et al. 1991; Tsurutani and Gonzalez. 
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1997; Richardson et al, 2002]. Nevertheless, recent stud- 
ies showed that some major storms may be driven by CIRs 
[Zhang et al. 2003; Richardson et al. , 2006], although their 
Dst values were not too far below -100 nT. Regardless of the 
solar origin, the geo-effective solar wind is usually a period of 
with prolonged and enhanced southward- directed magnetic 
field (B s ) that allows efficient solar wind energy transport 
into the Earth's magnetosphere [e.g., Dungey , 1961; Gonza- 
lez et al, 1994]. This enhanced B s field could be embedded 
within any part (front or rear) of ICMEs, SHs and CIRs [e.g., 
Grooker et al, 1992; Wu and Lepping , 2002; Huttunen and 
Koskinen , 2004; Richardson et al. . 2006]. 

Routine associations between ICMEs observed in geo- 
space and CMEs observed at the Sun became possible af- 
ter the launch of the SOHO spacecraft. Because of unfa- 
vorable launching directions and limited angular spans, the 
majority of CMEs do not intercept the Earth. However, a 
front-side halo CME, which appears as an expanding cir- 
cular feature surrounding the coronagraph occulting disk 
and thus likely has component moving toward the Earth 
along the Sun-Earth line, is likely to produce an ICME at 
the Earth [Howard et al. , 1982]. Comprehensive association 
work, based on a large number of CMEs and ICMEs con- 
tinuously observed over years, have been carried out [e.g., 
Lindsay et al , 1999; Gopalswamy et al, 2000; Cane and 
Richardson , 2003; Schwenn et al, 2005]. In general, based 
on existing solar and solar wind observations, one is able 
to make unique CME-ICME association for about half of 
all ICME events. However, reliable one-to-one associations 
for other ICMEs becomes more difficult, mainly because 
multiple activity at the Sun results in complex interplan- 
etary flows or compound streams [Gopalswamy et al, 2001; 
Burlaga et al, 2002; Zhang et al, 2003] or provides sev- 
eral plausible candidate associations. Further, a number of 
ICMEs, including those causing major geomagnetic storms, 
were found not to be associated with any identifiable front- 
side halo CMEs [Zhang et al, 2003; Schwenn et al, 2005] 

In this paper, our focus is to identify the solar and IP 
sources that lead to major geomagnetic storms. Our com- 
prehensive search for the sequence of events includes the so- 
lar surface sources, flare activity, CMEs, ICMEs and CIRs. 
Various tracking methods are used to address not only the 
obvious one-to-one events, but also to provide the possi- 
ble sequences for all complex events and problem events as 
well. While the evolution of an event is from the Sun to 
the Earth, it is practical to work backward from the Earth 
to the Sun for reliable identifications. The organization of 
the paper is as follows: Section 2 discusses the selection of 
the major geomagnetic storms. In section 3, we describe the 
methods used to identify the IP and solar sources of these 
geomagnetic storms. In section 4, we list the properties of 
the identified solar and IP sources and discuss the statistical 
results. Section 5 summarize the results of this paper. 
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2. Selection of Major Geomagnetic Storms 

The workshop focused on the major geomagnetic storms 
that occurred between January 1996 and December 2005. 
This 10-year period extends from the start to most of the 
descending phase of solar cycle 23, which had two sunspot 
maxima in 2000 and 2001. The Dst index is a measure of 
the strength of the ring current and widely used for measur- 
ing the intensity of geomagnetic storms. We defined a ma- 
jor geomagnetic storm as a minimum in hourly Dst index 
falling below -100 nT, A similar threshold for major /intense 
storms has been used by other authors [e.g., Tsurutani et al . , 
1997]. Other indices may be used, such as the Kp index [e.g., 
Gosling et al, 1991; Richardson et al. , 2002], Further, if a 
period of high activity showed multiple Dst < —100 nT min- 
ima, we arbitrarily assigned these to a single storm event if 
the minima were separated by less than 24 hours, rather 
than define each minimum as a separate storm (except the 
two storms that occurred at 12 UT, August 6 and 06 UT, 
August 7, 1998, which corresponded to two well separated 
ICMEs). As will be noted later, both single and multiple 
solar CMEs were found to be responsible for minima within 
a “single” storm event. 

We identified 88 major geomagnetic storms in total from 
January 1996 to December 2005, using the selection crite- 
ria described above. The events through 2003 are based 
on the final Dst index, whereas those in 2004 and 2005 
axe based on the provisional Dst index, so it is possible 
that they may be adjusted slightly based on the final in- 
dex. (Dst data are obtained at http://swdcdb.kugi.kyoto- 
u.ac.jp/dstdir/index.html). The 88 storms are listed in Ta- 
ble 1, where the first three columns indicate the event ref- 
erence number, the storm peak time and the minimum Dst 
value, respectively. The other columns, which will be ex- 
plained later, describe the parameters for the solar and IP 
sources. 

Figure 1 shows the distributions of the storm strength 
(top panel), yearly occurrence rate (middle panel) and oc- 
currence rate as a function of calendar month (bottom 
panel). A majority of these events (60 out of 88; 68%) 
have minimum Dst between -100 nT and -150 nT. A fur- 
ther 10 events (11%) have minimum Dst between -150 nT 
and -200 nT. There are 18 “severe” storms (21%) with 
minimum Dst < —200 nT. The largest geomagnetic storm 
(Dst = —422 nT) occurred on November 20, 2003 [Gopal- 
swamy et al, 2005]. The yearly major storm occurrence 
rate was highest (~ 13 events per year) during 2000 - 2002 
around the time of maximum sun spot number (SSN). The 
occurrence rate was lowest in 1996 at solar minimum. The 
bottom panel of Figure 1 shows that the occurrence of major 
storms in general follows the well known semiannual varia- 
tion of geomagnetic activity [e.g., Russell and McPherron, 
1973; Oliver et al, 2002], that is, higher activity during the 
equinoctial months and lower activity around the solstitial 
months. The number of major storms peaks in April-May 
and in October-November, and is lowest in June and in De- 
cember (when no storms occurred). Interestingly, the num- 
ber of major storms around the fall equinox is almost twice 
that at the spring equinox, and there are 55 events during 
the second half of year compared with only 33 during the 
first half. 

3. Methods of Identifying Solar-IP Sources 
of Major Storms 

3.1. Identifying and Characterizing the IP Sources 

The primary physical mechanism for energy transfer from 
the solar wind to the magnetosphere is magnetic reconnec- 
tion between the IMF and the Earth’s magnetic field. The 



efficiency of this process mainly depends on the strength of 
the southward IMF. or more accurately, the dawn-dusk (— y) 
component of the electric field (E s= —V x B) [e.g., Dungey , 
1961; Perreault and Akasofu, 1978; Tsurutani and Gonza- 
lez. 1997]. One formulation for the Dst index [0 } Brien and 
McPherron , 2000] relates the (pressure-corrected) Dst * in- 
dex to the solar wind driver given by VB S , where VB S is the 
rectified value of V B z that is positive when B z is southward 
and zero when B z is northward. The equations are: 



= Q(VB S ) - 


Dst* 

r(VB s y 


a) 


Q(VB S ) = 


a(VB s - E c ) VB S > E c , 
0 VB S < Ec, 


( 2 ) 


t(VB s ) = Too e xp(^~— g-)- (3) 

The rate of change of Dst* is assumed to be proportional 
to VB S ( Q representing injection into the ring current) less 
a loss term represented by the recovery time r that depends 
on the strength of the ring current and is assumed to be 
proportional to Dst*. 

Since storms are driven by the solar wind magnetic fields 
and plasma impinging on the Earth, we used in-situ solar 
wind plasma and magnetic field observations from the Ad- 
vanced Composition Explorer (ACE) and WIND spacecraft 
to identify the IP sources of the geomagnetic storms in this 
study. For ACE data, covering events during early 1998 - 
2005. 64-s resolution data were examined. We also exam- 
ined solar wind ion composition data from the ACE/SWICS 
instrument. ACE is in orbit at the upstream LI point, so 
there is typically a ~ 20 — 60 minute delay for solar wind 
structures to transit from ACE to the Earth. For WIND 
data, 92-s resolution data were used. During the period of 
this study, WIND spacecraft executed a complicated trajec- 
tory in the near- Earth solar wind with a variable solar wind 
transit time delay of typically less than 1 hour. Because of 
the near-complete observations provided by two spacecraft 
together, we were able to deduce the IP sources for all 88 
major geomagnetic storms studied. 

Based on their plasma and magnetic signatures, we 
identified various types of structures in the near-Eath 
solar wind in association with the geomagnetic storms. 
These include ICME-driven shocks, SHs, ICMEs in- 
cluding magnetic clouds (MCs)(a sub-set of ICMEs in 
which the magnetic field is enhanced and rotates through 
a large angle [Klein and Burlaga , 1982]), and CIRs. 
To assist in these identification, we referred to the 
several existing catalogs. For shocks, we used the 
WIND shock list compiled by J. Kasper (MIT) ( http : 
/ /space. mit.edu/home/jck/shockdb/ shockdh.html), and 

the ACE shock list compiled by C. W. Smith (UNH) ( http : 
/ / www— ssg.sr.unh.edu / mag / ace /ACElists / obsiist.html) . 
For ICMEs, we referred to an updated version of the “com- 
prehensive” ICME list compiled by Cane and Richardson 
[2003]. In addition, we used lists of MCs and “cloud- 
like” ICMEs compiled by R. P. Lepping and C.-C. Wu 
{http : / /lepm fi.gsfc.nasa.gov/mfi/MCLl.html) [Lep- 
ping et at. 2005] and the magnetic cloud list of Huttunen 
et at [2005]. Considering plasma composition and charge 
states, we used the list of high Fe-charge state intervals that 
are frequently associated with ICMEs, compiled by Lepri 
et at [2001], supplemented by information on compositional 
and charge state anomalies, also typically associated with 
ICMEs, based on the study of Richardson and Cane [2004]. 

The storm of July 27, 2004 (Event 75 in Table 1) serves 
to illustrate the method of source identification, as shown 
in Figure 2. The top panel shows the Dst index, indicating 


that this storm had a minimum value of Dst = —197 nT at 
14 UT. The following five panels show time profiles of the 
IMF strength and north-south (z) component, velocity, pro- 
ton density, proton temperature and calculated plasma /3, 
respectively. The three solar images at the bottom will be 
explained later. The IP driver of the main phase of the storm 
was evidently the extended interval of southward magnetic 
field reaching values of ~ 20 nT that started at ~ 05 UT 
on July 27, and lasted for about 10 hours. There was also a 
separate interval of southward field from ~ 22 UT on July 
26 to ~ 02 UT on July 27 that depressed Dst just below 
-100 nT at ^ 3 UT. Dst then recovered in response to a 
northward turning of the IMF; note the ~ 2 hour delay in 
the Dst response due to the solar wind transit time from 
ACE and magnet ospheric effects. 

Examining the broader context of the solar wind driver, 
we identified the passage of a fast forward IP shock at 
22:27 UT (at ACE; 22:25 UT at WIND) on July 26 (indi- 
cated by the vertical red line in Figure 2), characterized by 
abrupt jumps in the solar wind magnetic field, speed, den- 
sity and temperature. The shock was followed by a “sheath” 
of shocked IP plasma characterized by enhanced, fluctuating 
field strength, speed, density and temperature, extending for 
about 4-hours. 

The interval between the two blue vertical lines is the 
probable time of passage of the ICME that was driving this 
shock. The signatures of ICMEs have been discussed ex- 
tensively [e.g., Neugebauer and Goldstein , 1997; Wimmer- 
Schweingruber et at, 2006; Zurbuchen and Richardson. 
2006]. Here, we note the abnormally low proton temper- 
ature, depressed below the expected temperature for nor- 
mal solar wind [Richardson and Cane, 1995] overlaid in 
red, together with the enhanced magnetic field, smooth ro- 
tation in field direction (evident in B z ), and low plasma j3 
that is characteristic of a MC. Other signatures (not shown 
here) include enhanced oxygen charge states observed by 
ACE/SWICS and bidirectional suptrathermal electron flows 
observed by the ACE solar wind plasma instrument. Thus, 
the extended region of southward field driving the main 
phase of this storm was associated with the passage of a 
MC. The short period of southward field producing the ini- 
tial phase of the storm was associated with the sheath of 
shocked plasma ahead of the MC. Compressed magnetic 
fields in sheath regions may be draped over around the ap- 
proaching ICME [e.g. Gosling and McComas. 1987]. This 
may lead to strong out of the ecliptic fields, perhaps ac- 
counting for the initial phase of this storm. Two notable 
features of this event are the high solar wind speeds, reach- 
ing ~ 1000 km/s, in the SH and MC, and the overall low 
solar wind densities compared to average values. 

Considering CIRs, regions of compressed plasma formed 
by the interaction of high-speed streams from coronal holes 
with the preceding slower solar wind, these can be recog- 
nized by their characteristic variations in plasma parame- 
ters, including enhancements in the magnetic field strength, 
plasma density, temperature, and flow deflections lying at 
the leading edges of corotating high-speed streams [e.g., 
Forsyth and Marsch, 1999] and references therein. Exam- 
ples of major storms in our study driven by CIRs have been 
illustrated by Richardson et al. [2006], so a sample event will 
not be discussed in the present paper. For a recent review 
of CIRs and associated geomagnetic activity, see the papers 
in Tsurutani et al. [2006]. 

3.2. Identifying Solar Sources 

To identify the solar sources of the IP structures such as 
ICMEs that drive the major storms studied, we predom- 
inantly used observations from instruments on the SOHO 
spacecraft. CMEs near the Sun are observed by the LASCO 



C2 and C3 coronagraphs, which have fields of view of 2 - 
6 R s and 4-30 R S: respectively [ Brueckner et al. . 1995]. 
There were LASCO observations for 80 of the 88 major 
geomagnetic storms studied. The 8 events with LASCO 
data gaps occurred mostly in 1998 and 1999 when SOHO 
lost control for many months. To identify the surface 
features of CMEs in the source region, observations from 
SOHO's Extreme- Ultraviolet Imaging Telescope (EIT) [De- 
laboudiniere et al, 1995]. which images the Sun's corona over 
the full disk and up to 1,5 R s , were used, in particular those 
in the 195 A passband which is dominated by Fe XII emis- 
sion and sensitive to a plasma temperature of about 1.5 MK. 
In addition to referring to the LASCO CME catalog gener- 
ated by NASA and The Catholic University of America in 
cooperation with the Naval Research Laboratory [Yashiro 
et al, 2004] (http : / /cdaw.gs fc.nasa.gov/CMEJist/). we 
also carefully examined all the LASCO and EIT images in 
a suitable period prior to each storm to search for any erup- 
tion features that might not have been included in the cata- 
log, and to confirm the nature of the cataloged events. The 
Michelson Doppler Imager (MDI [ Scherrer et al. 1995] pro- 
vided photospheric magnetograms. 

In addition to SOHO observations, we used “tra- 
ditional” synoptic data, such as daily NO A A so- 
lar event reports, which include data on soft X-ray 
flaxes, filament eruptions and active regions (http : 
/ lwww.sec.noaa.gov/ ftpdir/ indices/). These data com- 
plement and reinforce the SOHO LASCO/EIT observations. 
We have also used X-ray coronal images made by the Yohkoh 
Soft X-ray Telescope (SXT) [Tsuneta et al , 1991] while it 
was available (Yohkoh was permanently lost in December 
2001) to search for possible eruption signatures. X-ray imag- 
ing observations made by the Soft X-ray Imager (SXI) on the 
GOES satellites [Hill et al. 2005] have also been used when 
available. For events from February 2003 onward, observa- 
tions from the Solar Mass Ejection Imager (SMEI) [Jack- 
son et al ., 2004; Webb et al. 2006] were used to help track 
CMEs to larger distances from the Sun than is possible with 
LASCO. 

The method of identifying the corresponding solar source 
of an existing ICME is straightforward, albeit ambiguous in 
many cases, that is to find a front-side halo (partial or full) 
CME at a reasonably earlier time, which depends on the 
transit time of the CME from the Sun to the Earth [e.g. 
Webb et al. 2000; Zhang et al. 2003]. The justification of 
this method is that there must be a cause-and-effect rela- 
tionship between solar and IP events, even though current 
observations only cover the near-Sun space, through remote- 
sensing, and the near-Earth space through in-situ sampling. 
However, the lack of imaging observations in the vast region 
between the Sun and the Earth through which CMEs can 
travel for days without direct tracking, contribute to the 
ambiguity of any such associations. 

Among the many CMEs observed at the Sun, halo CMEs, 
seen as an expanding circular bright feature fully surround- 
ing the coronagraph occulting disk (augular width 360°), 
are believed most likely to hit the Earth [e.g., Howard et al. 
1982]. The large angular width observed is attributed both 
to the projection effect and a large intrinsic width, indicat- 
ing the CME axis is likely directed along the Sun-Earth line, 
either toward the Earth if originating from the front-side of 
the Sun, or away from the Earth if originating from the back- 
side of the Sun. In addition to ’'full” halo CMEs. we also 
consider “partial halo” CMEs (apparent angular width > 
120°) in the solar source identification. To verify the surface 
source region of a CME, we mainly use EIT observations, 
which often manifest the CME origin with several erup- 
tive features, including large scale coronal dimming [e.g., 
Thompson et al, 1998] and post-eruption loop arcade (the 
counterpart of the more familiar post-flare loop arcade in 
H a ). These eruptive features are often associated with lo- 
calized coronal bright enings (the counterparts of flares in 
EUV wavelength). 


Considering the complexity in associating CMEs with 
ICMEs, we exploited an iterative process with multiple 
steps. First, we found all candidate front-side halo CMEs 
within a 120-hour-long search window before the arrival time 
of the ICME-driven shock (or other upstream disturbance 
if there was no fully-developed shock, or the ICME arrival 
if there was no upstream disturbance). The 120-hour-long 
search window corresponds to a 1 AU transit speed of 347 
km/s and is large enough to cover most possible CMEs 
sources except for extremely slow events. The large search 
window may produce several CME candidates, but further 
steps help to distinguish between likely and unlikely asso- 
ciations. The next step is to reduce the search window by 
estimating the CME transit time based on in-situ solar wind 
velocities at the location of shock arrival. Since fast CMEs 
tend to decelerate when moving through the slower solar 
wind, this method will give an upper estimate for the travel 
time. This method is not applicable to slow ICMEs because 
the corresponding, initially slow, CME may be accelerated 
by the ambient solar wind. In such cases, the full 120-hour 
window is used, and this may even be extended if the ICME 
of interest is extremely slow. The third step is that, for each 
remaining candidate CME in the search window, we consider 
whether the CME speed at the Sun is consistent with the 1 
AU transit speed implied by an association with the 1 AU 
shock/ICME, and with the in-situ solar wind speed. 

We recognize that the observed CME speed projected on 
the plane of the sky may not directly indicate the earthward- 
directed speed. Nevertheless, these speeds tend to be loosely 
correlated. Comparison with statistical studies of the rela- 
tionship between CME speeds and 1 AU transit times, e.g. 
[Cane et al. 2000; Gopalswamy et al, 2000; Zhang ei al. 
2003; Xie et al , 2004; Schwenn et al. 2005] can help to 
indicate whether a given CME-shock/ICME association is 
plausible or unlikely. We also take into consideration the 
solar source location implied by the CME/eruptive features. 
For example a central meridian source might be favored over 
a near-limb source, in particular if an ICME or magnetic 
cloud is involved in generating the storm. We should em- 
phasize that the CME-ICME associations were considered 
by the working group members both individually (often us- 
ing variations on the approach outlined above, and taking 
into account additional information, such as energetic parti- 
cle observations which may link solar events and interplan- 
etary shocks) and collectively, to arrive at a consensus. 

We will use the storm on July 27, 2004 (Figure 2) as an 
example to illustrate the identification process. The solar 
wind speed at shock arrival is ~ 900 km/s. If we simply 
assume that the CME- driven shock travels from the Sun at 
this constant speed, a travel time of ^ 46 hour is implied, 
suggesting (since this is a “fast” event at 1 AU) an CME 
event after 00 UT, July 25 as the source. Examining the 
LASCO CME catalog as well as the related images, there 
was only one halo CME in the search window, at 14:54 UT 
on July 25. This had a high projected speed (1333 km/s) 
which was consistent with the fast ICME seen at Earth al- 
lowing for some deceleration in the inner heliosphere. A 
direct association can also be demonstrated for this event 
using energetic particle observations which show an increase 
commencing at the time of the CME [Cane et al , 2006] that 
reaches peak intensity in the vicinity of the passage of the 
ICME-driven shock. This CME was associated with an Ml.l 
soft X-ray flare located at N04°H / 30°. The eruption at the 
surface was accompanied by a coronal dimming as shown in 
the running-difference EIT image (bottom middle panel of 
Figure 2). This CME/flare originated in NOAA AR 0652 
as indicated in the MDI magnetogram (bottom left panel of 
Figure 2). 



We should stress that it is not sufficient to use the time of 
the storm peak together with a plausible 1 AU transit time 
(e.g., based on the observed solar wind speeds at the time of 
the storm) to estimate the time of the solar source. Rather, 
it is important to examine and characterize the solar wind 
structures within which the geo-effective region is embeded, 
and then estimate the source timing. The effect of this dis- 
tinction is illustrated by the event in Figure 2: the peak of 
the storm is ^ 16 hours after the arrival of the shock and 
~ 12 hours after the arrival of the MC. These intervals are 
a significant fraction of the 1 AU transit times of the shock 
and ICME. Another point to note before leaving this event 
is that the two Dst minima in this storm result from two 
geo-effective regions, in the sheath and MC. associated with 
a single solar event. It therefor should not be assumed that 
multiple Dst minima within a storm interval indicate that 
multiple solar events are involved. 

3.3. Storms Involving Complex Solar Wind structures 
and Multiple CMEs 

We classify the solar-IP drivers of the major geomag- 
netic storms into three broad categories: S-type, M-type 
and C-type. S-type events are storms caused by single 
CMEs/ICMEs such as the July 24. 2004 storm described 
above. M-type are caused by multiple CMEs/ICMEs as dis- 
cussed in this section. The C-type axe for storms caused 
by CIRs [ Richardson et al, 2006]. For an M-type event, 
the storm is associated with complex solar wind structures 
that appear to involve multiple SHs and/or ICMEs. Two or 
more CMEs interact with each other in IP space, produc- 
ing such complex flows [Burlaga et a/., 2002; Zhang et al, 
2003; Wang et al , 2003]. Direct observations of the inter- 
action between two CMEs near the Sun have been reported 
[ Gopalswamy et al. 2001]. The M-type events are treated 
as a separate category from S-type because of the apparent 
differences in terms of the propagation/arrival of ICMEs. 
the resulting IP structure and the geo-effective components. 

One interesting variety of M-type events that we have 
noted is when a storm is generated by a faster ICME-driven 
shock propagating into the trailing edge of a slower ICME 
that originated in an earlier event at the Sun. An exam- 
ple is the storm of November 8, 1998 (Event 16, minimum 
Dst — —149 nT) shown in Figure 3. This storm was clearly 
generated by the region of southward magnetic field between 
21 UT, November 7 and 05 UT, November 8- The ACE 
plasma and field data show a weak shock at 07:36 UT on 
November 7 followed by a probable ICME commencing at 
~ 21 UT and indicated, for example, by the low proton 
temperature (black shading), enhanced magnetic field in- 
tensity, and enhancement in the solar wind O 7 /O 6 ratio. 
The southward magnetic field in this structure generated 
the onset of the storm, reaching levels of Dst ~ —100 nT. 
A second, stronger shock, propagating through the ICME 
passed ACE at 04:21 UT on November 8. The magnetic 
field in the ICME was starting to turn toward the eclip- 
tic at this time. However, the combination of the shock 
compression, which doubled the magnetic field strength and 
prevented the southward field strength from decaying, and 
the increase in solar wind speed, enhanced the y-component 
of the solar wind electric field, thereby strengthening storm 
activity and producing the peak of the storm. We suggest 
that ICME-associated plasma forms the post-shock sheath, 
at least to the end of the interval shown. Note that the 
field here turned northward, causing Dst to decline rapidly 
after the storm peak. We associate the shock on November 
8 with a 1119 km/s halo CME with a source at N22°W18° 
on November 5. Often in such situations, the source of the 
slower shock/ICME is less easily established. In the case of 
the shock on November 7, however, we suggest that a 523 


km/s halo CME at 07:54 UT on November 4 originating 
from a quiet-Sun region associatd with a quiescent filament 
is a likely candidate. We classify this storm as M-type be- 
cause, although the arrival of the November 8 shock is clearly 
associated with the peak of the storm, the presence of the 
southward fields in the preceding ICME is also required to 
generate the storm. 

Before leaving this event, it is worth commenting on the 
chance juxtaposition of the November 8 shock, Earth and 
preceding ICME that generated the storm peak. Had the 
timing been slightly different, the storm peak strength could 
have been substantially different. For example, had the 
shock been delayed relative to the ICME by as little as an 
hour or so, it would have encountered a region of north- 
ward field. Hence, the shock-ICME interaction would not 
have contributed to the storm. If the shock had arrived 
an hour or two earlier, it would have encountered stronger 
southward fields in the ICME, and an even more intense 
storm might have been generated. This clearly illustrates 
that while for S-type events involving one CME, there may 
be some hope in the future of predicting the geoeffectiveness 
using solar observations to infer the CME magnetic field 
structure, a similar prediction is far more problematical for 
M-type events. 

4. Results and Discussions 

4*1. Table of Solar and IP Sources 

Based on the methods described above, we have identi- 
fied the solar and IP sources of the 88 major geomagnetic 
storms during 1996 - 2005. The results are summarized 
in Table 1. Columns 1 to 3 give the properties of each 
geomagnetic storm, as discussed earlier. In column 4, we 
list the overall solar-IP source type (S, M and C). Columns 
5-10 describe the properties of solar source, and columns 
11-14 the properties of the IP sources contributing to the 
geomagnetic activities. Column 15 indicates a (somewhat 
subjective) confidence level for our identifications, given as 
1 to 3 in descending order of confidence. In the final col- 
umn, U F” indicates that there are additional comments in a 
footnote. In many cases these summarize critical comments 
on the proposed associations or alternative proposals, from 
working group members. 

Considering the properties of solar sources, the time in 
column 5 refers to the first appearance of CME in the 
LASCO C2 coronagraph, the CME velocity in column 6 
refers to the average velocity of CME in the LASCO C2/C3 
fields of view, and the angular width in column 7 is the 
apparent angular span of the CME in the plane of the sky 
measured in the C2 field of view. These values we generally 
obtained from the on-line LASCO CME catalog. However, 
in a few cases, they refer to previously unlisted CMEs that 
were identified from a re-examination of the LASCO images. 
Column 8 shows the magnitude of the soft X-ray flare as- 
sociated with the source CME. Column 9 shows the surface 
source region type, whether an “AR” (active region) fol- 
lowed by the four digit NOAA AR number, “QS” (quiet sun 
region), or “CH” (coronal hole). Note that quiet Sun regions 
here refer to any region on the surface of the Sun outside 
the traditional active regions and coronal holes. As for the 
sources of CMEs, they are often associated with erupting 
quiescent filaments or filament channels. 

Column 10 gives the heliographic coordinates of the sur- 
face source region. This generally corresponds to the H a 
flare location reported by NCAA SEC. When no H a flare 
location is reported, we used EIT images to measure the 
source coordinates, given by the location of the compact 
brightening, if observed, or the center of the dimming region 
if no brightening was observed. If the surface source regions 
of the CME candidates are unknown, because of the absence 
of any clear eruptive signatures on the disk in images from 



EIT and other instruments, this is indicated by “UNK” in 
columns 9 and 10. For those events with LASCO/EIT data 
gaps, the solar source could still be identified in some cases 
(events 6, 13, 14, and 21) because a major long-duration so- 
lar flare occurred at an appropriate time (based on consid- 
eration of transit times and in-situ solar wind speeds, solar 
particle events, etc) and location. For these events, the time 
in column 5 is the flare onset time, followed by “(F)” to em- 
phasize that this is not a CME time. Otherwise, “DG” in 
these columns indicates a gap in LASCO and/or EIT obser- 
vations and that it was not possible to identify a probable 
source using alternative observations. 

In the case of M-type events, there are multiple rows for 
each event listing the multiple CMEs that may contribute 
to the observed 1 AU solar wind structures. In each case, 
the first row indicates what we suggest is the “principal” 
solar driver. In the case of C-type events, the definitions 
of the parameters in the solar source columns are slightly 
different because of the different nature of the source. The 
time in column 5 indicates the central meridian transit time 
of the centroid of the associated coronal hole measured from 
EIT images. The time is followed by “(CH)” in order to 
emphasize that this does not refer to a CME source. The 
heliographic coordinate in column 10 indicates the latitude 
of coronal hole centroid when it crosses central meridian. 

Considering the properties of IP sources, column 11 char- 
acterizes the solar wind components that contribute to the 
storm, while columns 12, 13 and 14 show the time of the 
CME-driven shock (or disturbance, at either ACE, indicated 
by “A”, or WIND by “W”), and the start and ending times 
of the ICME. In column 11, We indicate in bold typeface 
the specific component (s) that contains the peak of the ge- 
omagnetic storm. Normal type indicates that the structure 
contributes to enhanced geomagnetic activity, but only to 
levels of > —100 nT. A plus sign indicates a simple suc- 
cession of components, while a dash indicates an “interac- 
tion” between the components. For example, for event 3, 
a sheath and magnetic cloud contribute to the geomagnetic 
activity. The sheath does not drive Dst to major storm lev- 
els, while the magnetic cloud includes the peak of the major 
storm. In contrast in event 5, though the same structures 
$xe present, the sheath drives the peak of the storm. For 
M-type events, it can be difficult to summarize in a com- 
pact way, or even to identify unambiguously, the various 
components present, but the nature of the specific compo- 
nents driving the storm is indicated, e.g., “SH(M)” means 
the presence of a sheath- like region that may include features 
(such as additional shocks) that suggest that more than one 
solar/ interplanetary event contributes. The situation where 
a shock is running into a preceding ICME or magnetic cloud, 
as discussed in relation to Figure 3 is indicated by PICME- 
SH or PMC-SH, respectively. 

When identifying the solar sources, we have found that 
storms generated by a single slow ICME present a major 
challenge. In particular, it is perplexing that, for about 10 
events (events 2,7, 28,31,34,36, 40, 58, 66, and 76 in Ta- 
ble 1), we were not able to find any halo CME candidates 
in the plausible search window with any apparent surface 
signatures (flare, filament eruption, dimming, brightening, 
or loop arcade), even though all the current solar disk ob- 
servations. including EIT and SXT, were available. Similar 
“problem events” have been reported ealier [Webb et al . , 
2000; Zhang et al. , 2003], and more recently, Schwenn et al. 
(2005) reported that about 20% of ICMEs observed at the 
Earth, regardless of the intensity of the resulting geomag- 
netic activity, were not preceded by an identifiable front-side 
halo CME (see also Cane and Richardson [2003]) Neverthe- 
less, similar to the finding of Zhang et al. (2003), we could 
always (except for event 40) identify a slow halo CME that 
occurred 4 to 6 days before the arrival of the corresponding 
ICME and was consistent with the inferred transit times 


from the CME and ICME speeds. Since there were no sur- 
face signatures, these slow solar CMEs would conventionally 
be regarded as backside halo CMEs. However, one possibil- 
ity, as suggested by Zhang et al.(2003), is that such slow 
CMEs could be from the front side but originate high in the 
corona, yielding little response in the low corona. In Table 1, 
we report such slow CMEs with unknown surface sources in 
columns 9 and 10. We emphasize that such reported CME 
sources are highly ambiguous, and caution must be exercised 
for further study of these events. 

To give an idea of the confidence of the identifications, we 
have assigned for each event a confidence level (indicated in 
column 15 of Table 1). Levels 1, 2 and 3 indicate, with 
decreasing level of confidence, the most unambiguous, plau- 
sible, and ambiguous/uncertain identifications, respectively. 
For only 46 (52%) of the storms would we regard our as- 
sociations as “level 1”'. These include most of the S-type 
and C-type events. The 27 (31%) level 2 storms include 
most M-type events and a few S-type. There are 15 (17%) 
events in level 3, including all the 10 problem events men- 
tioned above. In the following sub-sections, we summarize 
the properties of the solar and IP sources listed in Table 1 
and discuss their implications. 

4.2. On the Types of Overall Solar-IP Sources 

In Figure 4, we show the distribution of the three solar- 
IP source types for the 88 major geomagnetic storms dur- 
ing 1996 - 2005. The total numbers of S-type, M-type and 
C-type events are 53 (60%), 24 (27%), and 11 (13%), re- 
spectively. Hence, nearly two-thirds of these major storms 
were generated by single events at the Sun, and around an- 
other quarter involved multiple solar events. Considering 
S-type and M-type events together, we conclude that 77 

87%)* of the major storms in our study were driven by 
ICMEs (including the related upstream SH) and hence orig- 
inated from eruptive solar events, the remainder being asso- 
ciated with CIRs and hence with coronal holes. This result 
agrees with previous studies that have concluded that ma- 
jor geomagnetic storms are predominantly caused by ICMEs 
and their related structures [Gosling et aZ., 1991; Tsurutani 
et al. 1997; Richardson et al. 2001]. 

Nevertheless,, we also want to stress the non-trivial frac- 
tion (~ 13%) of these major geomagnetic storms that were 
driven by CIRs. Detailed analysis of these events (9 events 
from 1996 - 2004) has been reported by Richardson et al 
[2006]. This is a somewhat surprising result but is also a 
consequence of the -100 nT Dst storm threshold chosen for 
the workshop - the strongest CIR- associated storm had a 
Dst minimum of -131 nT so all these events would have 
been excluded had a slightly lower Dst threshold been cho- 
sen. Furthermore, we note that three of the 88 major storms 
were generated by the interaction of a CIR with an ICME. 
These were events 22 (October 22, 1999; Dst = —237 nT), 
event 58 (October 1, 2002; Dst = —176 nT), and event 76 
(August 30, 2004 ; Dst — —126 nT). These three events 
have been classified as S-type in the table because it is the 
presence of the ICME that is critical to the generation of 
the storm. 

The year-by-year distribution of event types is shown in 
Figure 5. In 1996, the year of solar minimum, there was a 
single major storm driven by a CIR. Otherwise, during the 
rise, maximum and declining phases of cycle 23, the major 
storms were predominantly driven by ICMEs with S-type 
dominating over M-type. C-type events were observed in 
1996 and 1998, were absent during 1999 - 2001 around so- 
lar maximum even though low latitude coronal holes and 
their associated streams were still typically present [Luh- 
mann et al. 2002], and reappeared in 2002 through 2005 
during the declining phase of the cycle. The asymmetry in 
the number (3 versus 8) of CIR.- generated storms between 



the rising and declining phase of the cycle, with more dur- 
ing the declining phase, has been noted in other studies [e.g., 
Richardson et al , 2001]. Nevertheless, most major storms 
were still driven by ICMEs during 2002 - 2005. 

For the 77 CME-driven storms events, around two thirds 
(53; 69%) were S-type and one third (24; 31%) M-type. The 
ratio of the numbers of S and M-type events does not show 
any clear solar cycle variation. Although we might expect 
M-type events to be more prominent at higher solar activity 
levels because of the higher CME rate, M-type events oc- 
curred throughout the solar rising, maximum and declining 
phases, except in 1997, when all 5 events were S-type, and 
S-type storms are still the most frequent type around solar 
maximum. The lack of a solar cycle dependence in the oc- 
currence of M-type events may be due to the fact that, for at 
least half of the 24 M-type storms, the responsible multiple 
CMEs originated from the same active region rather than 
from separate solar source regions. Such “super” active re- 
gions may appear at any phase of the solar cycle. 

4.3. On Geo-efFective Solar Wind Components 

Column 11 of Table 1 indicates that several configurations 
of IP structures gave rise to the major storms. For S-type 
events, the ICME and/or the upstream SH can contribute. 
We find that the storm peak was driven by the sheath in 12 
of these events (22%), by an ICME that is a magnetic cloud 
in 30 events (57%) and by a non-cloud ICME in 11 events 
(21%). Hence, a majority of major storms involving a single 
CME/ICME were driven to storm maximum by a magnetic 
cloud. For the M-type events, the IP drivers are typically 
more complex, and involve multiple structures. Neverthe- 
less, in most cases the storm driver can be characterized. In 
rare cases, such as event 10, a single driver among the vari- 
ous structures that pass the Earth (in this case a magnetic 
cloud) can be identified. A more common situation is that 
the storm peak is driven by a SH region or an ICME region 
that appears to include multiple components (indicated by 
SH(M) and ICME(M) respectively) that presumably reflect 
the complexity of the solar source. Multi-component SH re- 
gion drive 9 storms and multi-component ICME or MC re- 
gions another 6 storms. The situation illustrated in Figure 3 
in which a storm is caused by a shock propagating through 
a preceding ICME. drives the peak of 9 M-type storms, and 
hence is responsible for ~ 10% of all 88 major storms in this 
study. 

Considering the 53 S-type and 24 M-type CME-driven 
storms together, the geo-effective components axe MCs in 
33 events (43%), ICMEs without clear cloud signatures ac- 
count for another 14 events (18%), SH regions for 21 events 
(27%), and, as noted above, shocks propagating through 
preceding ICMEs/MCs in 9 events (12%). Hence, consis- 
tent with other studies, MCs form the most important class 
of IP drivers of major geomagnetic storms [ Wu and Leaping, 
2002; Huttunen et al , 2005]. This is despite the fact that 
only a minority of ICMEs at Earth, in particular around 
solar maximum, have magnetic cloud signatures [Richard- 
son and Cane . 2004]. The reason is that the magnetic fields 
associated with magnetic clouds can, if correctly oriented, 
provide the extended intervals of strong southward field that 
drive major storms, such as in Figure 2. Other ICMEs typi- 
cally have less organized, more irregular magnetic fields that 
may also be less enhanced, and hence non-cloud ICMEs are 
typically less geoeffective. Nevertheless, even if a magnetic 
cloud is present, it may not drive the peak of the storm if 
the cloud field orientation is not conducive for storm gen- 
eration. For example, in event 5, it is the sheath ahead 
of the magnetic cloud that drives the peak of the storm. 
More than half of the major storms are associated with other 
structures which have less organized magnetic structure, and 


hence in principle have less “predictable” geomagnetic con- 
sequences [Huttunen and Koskinen , 2004]. 

4.4. On Solar CMEs Associated with Major Geomagnetic 
Storms 

Except for the ~ 10% of events driven by CIRs, all the 
other major geomagnetic storms in our survey were caused 
by IP transients following solar CMEs. After excluding 
events that occurred during LASCO data gaps, we were 
able to identify 68 CMEs that were the likely solar sources 
of these storms, as given in Table 1. Apparently, these 
68 CMEs are the most effective in producing geomagnetic 
storms among thousands of CMEs observed during 1996 - 
2005. When summarizing the properties of these CMEs, 
only the presumed “principle” CME (shown as the first 
CME in the list of possible multiple sources in the event 
table) is included for M-type events. 

Considering the apparent angular size of these CMEs, 46 
(68%) were full halo CMEs, and 22 (32%) were partial halo 
CMEs. Clearly, partial halo CMEs should be considered 
when searching for the solar drivers of major geomagnetic 
storms. During the same period, LASCO observed 1187 halo 
CMEs of which 378 (32%) were full halos, and 809 (68%) 
were partial halos. Comparing with the number of similar 
CMEs that produced major storms, we estimate that about 
one out of eight full halo CMEs (or one quarter of front-side 
full halo CMEs, assuming that around half of halo CMEs 
originate on the backside of the Sun) will cause a major 
geomagnetic storm, and about 1 in 36 partial halo CMEs 
will do. If all LASCO CMEs, 10410 in total in the period 
of interest, are considered, on average only 1 out of ~ 150 
CMEs will cause a major storm. Since halo CMEs comprises 
only a small fraction of all CMEs observed, it is practical to 
use these relatively rare events to predict the interception of 
an ICME by the Earth, and hence the possible generation 
of a geomagnetic storm. However, these is certainly not a 
one-to-one association between halo CMEs and ICMEs at 
Earth. About 15% front-side halo CMEs may not inter- 
cept the Earth, and some 20% of ICMEs are not preceded 
by identifiable front-side halo CMEs [Schwenn et al , 2005] 
furthermore, when an ICME does intercept the Earth, the 
magnetic field configuration still has to be conducive for the 
generation of a major storm. The ICME rate at Earth [ Cane 
and Richardson , 2003], far exceeds the rate of major storms, 
for example by a factor of ^ 4 around solar maximum. 

In Figure 6, we display the speed distribution of the 68 
CMEs associated with major geomagnetic storms. Remark- 
ably, the distribution has a wide range from ~ 60 km/s to 
~2800 km/s with evidence of a peak at about 900 km/s. The 
average (median) speed of the 68 CMEs is 945 km/s (875 
km/s). A similar average speed (855 km/s) was obtained by 
Gopalswamy [2006] for a set of 55 geoeffective CMEs. For 
comparison, the average (median) speed of all 10410 CMEs 
in the study period is 472 km/s (410 km/s), and the av- 
erage (median) speed of all 1187 halo CMEs is 767 km/s 
(636 km/s). Hence, the major storm-associated CMEs are 
on average around twice as fast as the all-CME average, in 
agreement with recent results [Webb. 2002; Yashiro et ak, 
2004]. Forty-five (66%) of the 68 major storm-associated 
CMEs have speeds in the LASCO C2/C3 fields of view that 
exceed 600 km/s. These properties are consistent with the 
expectation that major geomagnetic storms are usually due 
to fast halo CMEs, 

Nevertheless, the relatively small (~ 200 km/s) difference 
between the average speeds for all halo CMEs and major 
storm- associated CMEs suggests that strongly geo-effective 
halo CMEs cannot necessarily be distinguished from other 
halo CMEs on the basis of their speed alone, as discussed 
earlier by Zhang et al [2003]. Further, some very slow 
CMEs, though a small faction, can also generate major 
storms. Twelve (18%) of the 68 storm- associated CMEs had 



speeds of less than 300 km/s. These results emphasize the 
fact that speed alone is not the major factor determining 
geoeffectiveness. Rather, the configuration of the embed- 
ded magnetic fields is also important, as exemplified by the 
fact that most of these storms resulted from slow magnetic 
clouds at the Earth with speeds comparable to the ambient 
solar wind. 

Considering the association of major storms with GOES 
soft X-ray flares, we find that among the 77 CME- driven 
storms, 19 (25%) were associated with a X-class flare, 17 
(22%) with a M-class flare, 19 (25%) with a C-class flare, 
and 22 (28%) with either minor (B or A-class), or with 
no evidence of a flare. We conclude that, major (M or 
X-class) flares were associated with about one half of our 
major storms, and that around a third of the storms were 
not accompanied by a flare or only by a minor flare. There- 
fore using flares, the traditional indicator of solar activity, 
to predict geomagnetic storms is often far from satisfac- 
tory [ Gosling , 1993]. 

4.5. On the Solar Surface Source Regions Associated 
with Major Geomagnetic Storms 

Figure 7 summarizes nature of the solar surface source re- 
gions where the major storms in our study originated (Col- 
umn 9 of Table 1). For 3 of the 88 events, there were in- 
sufficient data (e.g,, data gap in LASCO/EIT observations, 
and no major flares reported in a plausible time window) for 
the source to be inferred. In the case of M-type events, we 
only include the source of the principle CME. We find that 
54 storms (~ 63%) originated in active regions, 10 (12%) 
in quiet Sun regions, and 11 (13%) were associated with 
coronal holes. Here, quiet Sun region is a general reference 
to any coronal region other than active regions or coronal 
holes. It should be noted though that even when a CME 
originates outside an active region, it is usually associated 
with a quiescent filament or filament channel overlying a 
magnetic inversion line in the photosphere. For the remain- 
ing 10 (12%) events we were unable to identify any solar 
surface signature and hence the nature of the source region 
is unknown. 

Thus, while half of the major geomagnetic storms orig- 
inated in active regions, a similar number originated out- 
side active regions. Nevertheless, active regions remain the 
source of the largest storms. The ten largest storms (mini- 
mum Dst < -271 nT) during 1996 - 2005 were all associated 
with active regions. For comparison, the largest storm that 
originated from a quiet Sun region reached Dst = —237 nT. 
Furthermore, the largest storm with an unknown surface 
source attained Dst = — 182 nT, and the largest storm from 
a coronal hole source had a minimum Dst of only -131 nT. 

In Figure 8, we show the heliographic distribution of the 
source regions (Column 10 of Table 1). This distribution 
includes the 64 CMEs with identified surface sources. The 
other 24 events are excluded because they were associated 
with coronal holes (11 events), or unidentified sources (10 
events), or occurred within solar data gaps (3 events). The 
source locations lie within 35° N to 58° S latitude both for ac- 
tive region (red symbols) and quiet Sun sources (blue sym- 
bols), and 61 of the 64 source regions (95%) lie within 30° 
from central meridian. A possible explanation is that CMEs 
originating from higher latitudes propagate into the high 
latitude region of the heliosphere and do not intercept the 
Earth. 

Considering the longitudinal distribution, 56 of the 64 
source regions (88%) lie within 45° from central meridian, 
(77%) within 30°, and 34 (53%) within 15°. Hence, the 
vast majority of major storms arise from solar sources that 
are close to central meridian. Nevertheless, the sources also 
show an east- west asymmetry that favors the western hemi- 
sphere and reinforces the similar result from the study of 
Zhang et ah [2003]. Specifically, the sources extend to 85° W, 


but only to 58° E, and 42 lie on the western hemisphere, 
compared with 21 on the eastern hemisphere (one event is 
at central meridian). Hence, the ratio of number of western 
to eastern sources is ~2:1. The average (median) longitude 
of all the 64 events studied is 11°W (8°W). Geo-effective 
CMEs could be from far western regions but not from far 
eastern regions. This east-west asymmetry seems to be a 
general feature of the ICMEs that intercept the Earth, re- 
gardless the strength of geo-activity [Wang et ah. 2002; Cane 
and Richardson , 2003]. One possible explanation is that this 
asymmetry results from the deflection of CME trajectories 
by the spiral IP magnetic field [Wang et ah , 2004], 

4.6. Implication for Forecasting Major Geomagnetic 
Storms 

What are the implications of this study for forecasting 
major geomagnetic storms using solar observations? First, 
there may be a misconception that a major geomagnetic 
storm must be caused by an unusually fast halo CME from a 
strong active region accompanied by various energetic erup- 
tive signatures (e.g., major solar flares). Except for a few 
super storms, this was not the case for many of the major 
storms (Dst < -100 nT). In fact, some of these storms were 
caused by moderate speed CMEs that may originate outside 
of active regions, as well as by CIRs associated with coro- 
nal holes as described earlier. A central reason is that the 
driving electric field y-component depends on both the solar 
wind speed and B s , but the dynamic range of B s is greater 
than that of the solar wind speed. However, a combination 
of CME speed and magnetic field in ICMEs seem to have a 
high correlation with Dst index [ Gopalswamy , 2006]. Fur- 
thermore, activity is suppressed when the IMF is northward, 
so a fast ICME with a predominantly strong northward field 
will not generate a major storm. The size of a storm also 
depends on the time variation of the southward field compo- 
nent. Thus, a relatively slow moving MC with an extended 
region of enhanced southward field (such as event 15) can 
generate a major storm. Hence, the speed of a halo CME 
alone is an inconsistent predictor of a major geomagnetic 
storm. A major advance would be able to “predict” the 
interplanetary magnetic field configuration at 1 AU, in par- 
ticular for S-type storms involving only one CME/ICME, 
based on solar observations, but this is difficult at present. 

In the case of storms that involve more than one 
CME/ICME, a complicating factor for forecasting is that 
it is the details of the magnetic structures formed by the 
interaction of these transients (and their associated shocks) , 
both with each other and with the ambient solar wind, that 
determine the resulting level of geomagnetic activity. The 
precise path of the Earth through the structure is also a fac- 
tor. Thus, it is unlikely that even a relatively complete MHD 
simulation of two CMEs launched towards the Earth would 
ever include sufficient information to be able to model the 
resulting fields at 1 AU on the necessary few-hour timescales. 
Information from upstream spacecraft would help to assess 
the likely geomagnetic impact, but the interacting structures 
may still evolve before reaching Earth. 

5. Summary 

We have investigated the solar and IP sources of the 88 
major geomagnetic storms ( Dst < —100 nT) that occurred 
during 1996 - 2005 with the aim of providing a reliable as 
possible list of associations that is intended to provide a basis 
of future studies by the LWS CDAW participants and oth- 
ers. By combining remote-sensing solar observations, in-situ 
neax-Earth solar wind observations, and the wide range of 
experience of the Working Group members, we were able to 



identify with reasonable confidence the chain of sources for 
about 83% (73) of these events, although the detailed one- 
to-one association could not be established for those com- 
plex events involving multiple CMEs and ICMEs. We are 
uncertain of the origin of the other 17% (15) of the storms, 
mainly because their driving CMEs were not associated with 
noticeable eruption signatures at the solar surface. Detailed 
parameters of the solar and IP sources for each of the 88 
major geomagnetic storms have been provided. The main 
results are as follows: 

(1) Based on the overall solar and IP properties, the 
sources can be divided into three broad categories: S-type, 
driven by single CMEs and thier IP counterparts; M-type, 
associated with multiple CMEs/ICMEs, and C-type due to 
CIRs driven by high speed streams from coronal holes. The 
total numbers of S-type, M-type and C-type events are 53 
(60%), 24 (27%), and 11 (13%), respectively. 

(2) Of the 68 LASCO CMEs associated with these ma- 
jor storms, 46 (68%) were full halo CMEs, and 22 (32%) 
were partial halo CMEs. Their speed have a wide range (60 
km/s to 2800 km/s). The average speed (945 km/s) is about 
twice as fast as the average for all LASCO CMEs. About 
half (47%) of these storm-associated CMEs were accompa- 
nied by major (X and M-class) flares. 

(3) For the 85 storms for which we could identify the solar 
surface source, we find that 54 (~ 63%) originated in active 
regions, 10 (12%) in quiet Sun regions associated with qui- 
escent filaments, and 11 (13%) were associated with coronal 
holes. The other 10 (12%) events originated from unknown 
surface source regions. 

(4) Major geomagnetic storms predominantly originate 
from sources near central meridian (e.g., 88% from with 
45° , and 77% from with 30° of central meridian) but show 
an east-west asymmetry with around twice as many storms 
originating on the western hemisphere than on the eastern 
hemisphere. 
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Figure 1 . Distributions of the minimum Dst (upper 
panel, bin size = 50 nT), yearly occurrence rate (middle 
panel), and occurrence rate per calender month (lower 
panel) for 88 major geomagnetic storms during 1996 - 
2005. The black curve overlaid in the middle panel shows 
the 1 80-day-running-average daily sunspot numbers in 
arbitrary units 





Stort Time (26-Jul~04 12:00:00) 



Figure 2. Geomagnetic, interplanetary and solar data 
related to the major geomagnetic storm (minimum Dst = 
—197 nT) on July 27, 2004 (event 75). The upper six 
panels, from top to bottom, show temporal profiles of 
the Dst index, solar wind magnetic field intensity (black) 
with the B z component (red) overlaid, solar wind veloc- 
ity, density, and proton temperature (black) overlaid with 
the expected temperature (red) (Richardson and Cane, 
1995), and the plasma (3. The solar wind data are from 
ACE in GSE coordinates. The solid and dotted blue ver- 
tical lines indicate the starting and ending times of the 
ICME, which in this case is a magnetic cloud. The verti- 
cal red line indicates the arrival time of the ICME-driven 
shock. The three images at the bottom, from left to 
right, indicate the source active region in a SOHO/MDI 
magnetogram for July 25, the coronal brightening accom- 
panying the associated CME observed by EIT (running 
difference image), and this CME shown in a LASCO C2 
coronagraph running difference image. 
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Figure 3. Geomagnetic and interplanetary data for the 
major geomagnetic storm (minimum Dst = — 149 nT) on 
November 8, 1998 (event 16). The panels show, from top 
to bottom, the observed Dst (black) with the predicted 
Dst index using the O’Brien & McPherron formula [2000] 
overlaid in red, the magnetic field intensity (black) with 
B z overlaid in red, the Y-component of the solar wind 
electric field, the solar wind velocity, density, proton tem- 
perature (black) and expected proton temperature (red) 
with the shaded black shading indicating where the pro- 
ton temperature falls below the expected temperature, 
Helium/proton ratio, and O 7 /O 6 ratio. The two verti- 
cal green lines indicate the arrival times of ICME-driven 
shocks. Here, the peak of the storm is caused by an in- 
terplanetary shock (~ 04 UT on November 8) propagat- 
ing through a preceding ICME which has an embedded 
strong southward magnetic field. 
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Figure 4. Distribution of the three types of solax-IP 
sources for the 88 major geomagnetic storms during in 
1996 - 2005. 





Figure 5. Solar cycle variation of the occurrence rate 
of the three types of solar-IP sources for the 88 major 
geomagnetic storms during 1996 - 2005. 




Figure 6. Distribution of the plane of the sky speeds for 
the 68 CMEs observed by SOHO/LASCO that resulted 
in major storms. 
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Figure 7. Types of solar surface source regions for the 88 major geomagnetic storms during 1996 - 2005. 




Surface Source Regions of Major Geomagnetic Storms 



Figure 8. Heliographic locations of the 64 identified sur- 
face source regions for the CMEs that resulted in major 
geomagnetic storms during 1996 - 2005. 




